In-situ hydrothermal synthesis of CNT decorated by nano ZnS/CuO for simultaneous removal of acid food dyes from binary water samples

The zinc sulfide/copper oxide–carbon nanotube nanocomposite (ZnS/CuO-CNT) was fabricated by using an in-situ hydrothermal synthesis method and was used for simultaneous ultrasound-assisted adsorptive removal of a binary mixture of ponceau 4R (P4R) and tartrazine (TA) acid food dyes from contaminated water. The as-synthesized ZnS/CuO-CNT was described by FESEM, XRD, FTIR, BET, and zeta potential analysis. The results included nested network morphology, high purity with the crystalline structure, oxygen-containing functional groups, mesoporous/micropores texture with cumulate interspace, specific surface area of 106.54 m2 g-1, and zero-point charge (pHzpc) of 5.3. In adsorption experiments, the simultaneous effect of main independent variables, including solution pH, adsorbent dosage, concentration of each dye, temperature, and sonication time on the removal efficiency of dyes was studied systematically using the central composite design (CCD) method based on response surface methodology (RSM). Also, the second-order multivariate equation was presented to determine the relationship between the removal efficiencies of P4R and AT dyes and six independent effective variables. The high correlation coefficient (R2 ≥ 0.99), significant p-value (P < 0.0001), and non-significant lack-of-fit (P > 0.05) showed the high accuracy, and validity of the proposed model to predict the removal efficiency of P4R and TA acid food dyes. The experimental removal efficiency for P4R and TA dyes was found to be 98.45 ± 2.54, and 99.21 ± 2.23, respectively. Also, the Langmuir maximum adsorption capacity for P4R and TA dyes was determined to be 190.1 mg g-1 and 183.5 mg g-1, respectively. Finally, the adsorbent's reusability was tested for six periods and could be reused repeatedly without significant reduction in adsorption performance.

respiratory illnesses 9 . Research has also shown that daily consumption of foods containing food additives has also been found to lower children's IQ by five points, according to research 10 . Ponceau 4R (P4R) and Tartrazine (TA) are two anionic synthetic food dyes that are always used together to form an orange dye 11 . P4R and TA have been widely used to make groceries more attractive and tasty because of their exceptional benefits such as low cost, high solubility in water, down microbiological pollution, and high consistency 12 . Nonetheless, P4R and TA belong to the azo dyes (-N = N-) with an aromatic and carcinogenic structure that can lead to side effects such as allergies, cancer, and neurological behaviors if consumed in excess [13][14][15] . In scientific references, different removal methods such as membrane methods 16 , electrochemical and sonochemical 17 , ion exchange resins 18 , coagulation 19 , flocculation, 20 nanofiltration 21 , and reverse osmosis 22 have been used to remove synthetic dyes. However, these methods have limitations such as high cost, sludge formation, or by-products limiting their use 23 . With the advancement of science and technology, the adsorption process has recently been regarded for the treatment of different organic compounds as well as heavy metal ions in contaminated water owing to its excellent qualities such as low raw material costs, ease of design, safety, and efficiency 24 . Researchers have studied the removal of dyes with different adsorbents 25,26 . Among the adsorbents used, carbon nanotubes (CNTs) are widely used to remove various contaminants due to their high surface area to volume and high adsorption capacity 27 . Also, different types of metal oxides/sulfides nanoparticles have been combined with carbon nanotubes to remove a variety of contaminants from wastewaters 28,29 . Recently, applications of ultrasound for wastewater treatment have been conducted by many researchers in ultrasound-assisted adsorption 30 and can be used as a very effective technique as an accelerating agent in increasing the adsorption capacity. The increased adsorption capacity in the presence of ultrasound is defined due to the effect of ultrasound in reducing the agglomeration of nanoparticles, increasing kinetics, and enhancing the affinity between adsorbent and adsorbate 31,32 . Nowadays, many statistical experimental design methods are applied for the optimization of process parameters. In the conventional one-variable-at-a-time statistical method, the optimization of a multivariable system follows one factor at a time and does not represent the combined effect, and is not able to study the interactions and complexities between independent variables 33 . On the other hand, the conventional optimization technique requires more experiments to determine the optimum conditions and takes prolonged time 34 . However, the experimental design method based on the central composite design (CCD) method along with response surface methodology (RSM) is a mathematical and statistical technique that is able to optimize and model the complex relationships between variables and system responses with a limited number of experiments 35 .
In this study, the carbon nanotubes were decorated by nano zinc sulfide/copper oxide (ZnS/CuO) through the in-situ hydrothermal method to fabricate the zinc sulfide/copper oxide-carbon nanotube nanocomposite (ZnS/CuO-CNT) and was used as an efficient adsorbent for simultaneous ultrasound-assisted removal of dye mixture P4R and TA from contaminated water. The structural characteristics of adsorbent was studied in detail by field emission scanning electron microscopy (FESEM), X-ray diffraction (XRD), Fourier transform infrared spectroscopy (FTIR), Brunauer-Emmett-Teller (BET), and zeta potential (ZP) analysis. The adsorptive removal process was performed on a laboratory scale based and optimized by the central composite design (CCD) method based on response surface methodology (RSM). The effects of independent variables such as solution pH, adsorbent dosage, initial concentration of each dye, sonication time, and temperature at different levels were evaluated. Finally, the adsorption isotherms, kinetics, thermodynamics, reusability, and efficiency in real samples were investigated.
Preparation. Functionalization of carbon nanotubes. The chemical oxidation route was carried out by HNO 3 acid treatment as follows 36 : 100 mL HNO 3 (15 mol L -1 ) was prepared and added into a round bottom flask. Then, 100 mg of pristine carbon nanotubes were dispersed into the acid solution by ultrasonication for 30 min and then refluxed for 5 h in an oil bath at 120 ºC. After centrifugation at 5000 rpm for 10 min, the HNO 3 -treated CNTs were recovered and washed with ultrapure water until the pH of the CNT solution neared neutral. Finally, the functionalized carbon nanotubes (CNTs) were dried in the oven at 120 ºC overnight.
Synthesis of CuO-CNT nanocomposite. The CuO-CNT nanocomposite was synthesized by the in-situ hydrothermal method as follows: The 0.7 g of dried CNTs were suspended in 50 mL of ultrapure water by sonicating  (Table 1) and in random order. In all experiments, 50 ml of a mixture of dyes P4R and TA was used at room temperature. For quantitative measurements through spectroscopy, the maximum wavelength of species adsorption in solution must be determined. Therefore, the zero-order absorption spectra of the solution of P4R and TA dyes, each with a concentration of 20 mg L -1 and a binary mixture of dyes in the range of 250 to 750 nm, were determined ( Fig. 1). Because the absorption spectra of P4R and TA dyes overlapped and displayed interference between their zero-order spectra, direct absorbance measurements could not be used to quantify their amounts simultaneously. First-order derivative spectrophotometry was created to tackle this challenge and determine the concentration of P4R and TA in binary mixes (Fig. 2). The use of zero-order derivative spectroscopy to separate overlapping signals and decrease the influence of spectral interposition caused by the attendance of another ingredient in a sample could be beneficial 38 . By detecting the absorbance signal at the first-order derivative wavelength, the residual concentrations of P4R and TA dyes in the binary mixture were monitored aided the simultaneous removal of dyes. The method of zero-crossing was used to choose two wavelengths for two dyes. The proper wavelength of P4R is 552 nm, for which the first-order derivative of TA is zero.
According to the first-order derivative method, the correct wavelength of TA is 374 nm (Fig. 2). As a result, the calibration curve for P4R and TA dyes was determined using the dye's first-order derivative spectrum at specified wavelengths (Fig. 3). The absorbance of the solution was then determined by UV-Vis spectroscopy at the specified wavelengths and the removal percentage was calculated by the below equation: In this equation, C 0 is the concentration of dye in the initial solution and C t is the concentration of dye in the final solution.
The accuracy of the first-order derivative method was determined by evaluating the recovery percentage by considering a binary mix of P4R and TA dyes with different concentrations. Table 2 shows a summary of the results. The recoveries and errors for the designation of the P4R and TA dyes in the binary mix using the first-order derivative spectroscopic technique were 96-103 percent and − 3.17 and + 3.62 percent, respectively, as shown in Table 2. Based on recovery perusal, the first-order derivative spectroscopic technique can correctly identify the amounts of P4R and TA dyes in the binary mixture 39 .

Experiment design.
To evaluate the performance of the ZnS/CuO-CNT adsorbent, a discontinuous system and experimental design method were used. For this purpose, the CCD-RSM design was applied for five independent variables with five levels. The solution pH, adsorbent dosage, initial concentration of each dye, ultrasound time, and temperature were considered the independent variables and removal percentage was considered dependent variables. Table 1 shows the independent variables, their range and experimental levels. In general, in this study, the number of experiments was 32. The results of experimental work in different periods were used to prepare a predictable model. This equation is a polynomial regression model discussed in previous studies 40 .

Results and discussions
Characterization. The surface morphology of the samples was investigated by the FESEM images and are shown in Fig. 4. In the FESEM image of CNT (Fig. 4a), nanotubes are seen as filaments of micrometer length and nanometer thickness that are densely intertwined with an almost smooth surface. In the FESEM image of CuO-CNT (Fig. 4b), it was observed that the CuO nanoparticles with a near-spherical shape and average diameter of about 20-40 nm were formed in the CNTs network. Also about ZnS/CuO-CNT (Fig. 4c), the growth of almost www.nature.com/scientificreports/ spherical nanoparticles with a size of less than 50 nm in the between and surface of the nanotubes can still be seen, which has created a nested network morphology. The crystal phases and purity of samples were analyzed by XRD patterns and are shown in  42 . The sharp diffraction peaks indicated that the CuO nanoparticles were grown with a high crystalline structure in the presence of nanotubes. The XRD pattern of ZnS/CuO-CNT (Fig. 5c) showed the appearance of new diffraction peaks compared to the XRD pattern of CuO-CNT, which attributed to the (111), (220), and (311) crystal planes of cubic ZnS phase (JCPDS no. 80-0020) 43 . Also, no additional diffraction peaks were seen indicating that the synthesis of ZnS/CuO-CNT was successful. The FTIR measurements were conducted to evaluate the surface functional groups of samples and the results are shown in Fig. 6. In the CNT spectrum (Fig. 6a), the bands at 610 cm -1 , 1100 cm −1 , 1565 cm -1 , and 1640 cm -1 can be assigned to C-H (out of plane bending), C-O, C = C (aromatic rings), www.nature.com/scientificreports/ and C = O (-COOH) stretching vibrations, respectively 44 . In the FTIR spectrum of CuO-CNT (Fig. 6b), the presence of Cu-O absorption bands along with CNT bands was verified at wavenumber 480-580 cm −145 . In the FTIR spectrum of ZnS/CuO-CNT (Fig. 6c), the Zn-S stretching vibrations were assigned at 497 and 998 cm −1 , which confirmed the successful formation of ZnS in presence of CuO-CNT composite 46 . Also, the absorption peaks related to absorbed water molecules around 1640 cm -1 were not shown in the figure due to the possibility of overlapping with C = O absorption peaks. The textural properties of the ZnS/CuO-CNT were studied by nitrogen adsorption-desorption isotherms (Fig. 7a) and defined as a type II isotherm with H3 hysteresis  www.nature.com/scientificreports/ loop, which indicates the presence of both mesopores and micropores texture with cumulate interspaces in the sample 47 . The BET model also estimated the specific surface area, total pore volume, and mean pore diameter of ZnS/CuO-CNT, which were 106.54 m 2 g -1 , 0.283 cm 3 g -1 , and 10.67 nm, respectively. The zeta potential of ZnS/CuO-CNT was determined as a function of pH, and the result are presented in Fig. 7b. It is observed that the zeta potential of ZnS/CuO-CNT in the range of pH 4.0 to 8.0 is more strongly dependent on pH and its value becomes from + 12 to − 32 mV as the solution becomes more alkaline. Also, the zero value of the zeta potential (isoelectric point) was obtained and the pH of the zero point charge (pH zpc ) of the ZnS/CuO-CNT was 5.3.
Modeling and optimization of adsorption. The modeling and optimization of adsorption process was performed by the CCD-RSM design, and adequacy of the model was tested through the analysis of variance (ANOVA). The ANOVA results based on the quadratic model are presented in Table 3.
The coefficient R 2 represents the ratio of the total changes in the response predicted by the model, which represents the ratio of the sum of squares due to regression (SSR) to the sum of squares (SST). A large R 2 and close to 1 are desirable, and a favorable agreement with the essential R 2 -Adj. The large R 2 confirms the satisfactory  www.nature.com/scientificreports/ conformity of the test data to the quadratic model 48 . F-value is an indicator to evaluate the significance of the model that the larger the numerical value, and the, more significant the model and P-value are less than 0.05 approved by the model 49 . The value of P-related factors greater than 0.1 do not have a significant effect on the model and by removing them, the model can be modified. As shown in Table 3, the fitting weakness test for the quadratic polynomial model fitted to the surface response data is not significant, indicating a successful fit of the data. Therefore, high correlation coefficient (R 2 ≥ 0.99), significant p-value (P < 0.0001), and non-significant lack-of-fit (P > 0.05) indicate high accuracy, and validity of the proposed model to predict the efficiency of simultaneous removal of dyes mixture by ZnS/CuO-CNT from aqueous solutions. Multi-nominal regression modeling between the response variable (removal %) and the corresponding coded values (X 1 , X 2 , X 3 , X 4 , X 5 , and X 6 ) of five different variables, respectively (solution pH, ultrasound time, adsorbent www.nature.com/scientificreports/ dosage, temperature, initial concentration of each dye) was created and finally, the best-fitted model was obtained as follows: Negative coefficients imply a diminishing influence of factors on dye adsorption, whereas positive coefficients suggest an increasing effect of variables 50 .
Also, the actual values plotted against the predicted values derived from the CCD model of dyes removal efficiency are shown in Fig. 7c. As can be seen, the plot indicated a good agreement between the actual and predicted values for the removal efficiency of dyes, which indicated the efficiency and accuracy of the predicted model. In   www.nature.com/scientificreports/ addition, the residuals need to be observed to ensure that the ANOVA assumptions are valid. The residual error analysis was studied by the residuals vs observation order plot as shown in Fig. 7d. It shows that the residuals in the plot fluctuate around the centerline and don't explicitly show a trend. Therefore, the randomization in the experiment was vindicated suggesting that the ANOVA assumptions are valid.
Simultaneous effect of independent variables on dyes removal efficiency. The simultaneous effect of independent variables such as pH (2.0-10), adsorbent dosage (5-15 mg), initial concentration of each dye (10-50 mg L -1 ), sonication time (2-16 min), and the temperature (5-45 °C) on P4R/TA dyes removal efficiency was displayed by the three-dimensional response surface diagrams as shown in Fig. 8. As can be seen, three-dimensional response surface diagrams were performed as a function of two independent variables that change in laboratory ranges and the third variable was the removal efficiency of the dyes, which was reported as a removal percentage (R%). The simultaneous effects of sonication time and solution pH on R% of P4R and TA are shown in Fig. 8a,b, respectively. It was observed that the R% has increased at lower pH ranges as well as longer sonication times in both dyes. The P4R and TA dyes are anionic due to the presence of sulphonic groups in their structures. Also, in acidic conditions, the adsorbent surface has a positive charge due to its zeta potential and pH zpc . Therefore, the generated electrostatic attraction can play an effective role in the adsorption process and increase the removal efficiency in acidic conditions. In contrast, with the increase of pH and more became alkaline, the positively charged sites of the adsorbent get decreased and trends to create the negatively charged. This situation leads to electrostatic repulsion between the negatively charged surface of the adsorbent and anionic dye causing a decrease in adsorption of the dyes from the aqueous solution. Similar findings were also reported in other studies 51,52 . Further, ultrasonic waves contribute to faster mass transfer processes, increase the diffusion coefficient, and also increase the chemical reaction of particles in solution by increasing the effective collisions interparticle collisions 53,54 . Therefore, the removal efficiency was increased with increasing the sonication time.
On the other hand, the equilibrium adsorption time was very fast in occurred in 12 min. Rapid adsorption is due to the surface bonding of the species to the adsorbent surface. It is believed that when the adsorption process is rapid, the speed limiting step is related to the transfer process in the liquid phase, such as penetration into the www.nature.com/scientificreports/ solution mass, penetration from the vicinity of the adsorb film to solid particles, or penetration into cavities filled with solution 55 . The simultaneous effects of sonication time and adsorbent dosage on R% of P4R is shown in Fig. 8c. As can be seen, with increasing sonication time and adsorbent dosage the removal efficiency of dye has increased until it has finally reached a state of equilibrium. The rate of rapid adsorption in the first minutes of the reaction can be interpreted as due to the presence of very high adsorption sites and the high concentration gradient of contaminant molecules. Gradually, however, the adsorption sites become saturated and the adsorption rate is controlled by the phenomenon of transfer from the outside to the inside of the adsorbent. Also, it is obvious that with increasing the amount of adsorbent, the number of available vacant active sites on the adsorbent increases compared to dye molecules and more adsorption sites remain unsaturated during the adsorption process causing the removal efficiency was increased 56 . In the adsorption process by ZnS/CuO-CNT, the simultaneous effect of adsorbent dosage and initial dye concentration on R%TA is illustrated in Fig. 8d. It has been shown that the R%TA increased with decreasing TA dye concentration and increasing the adsorbent dosage. In fact, in high dye concentrations, the vacant sites of the adsorbent are saturated and increase competition in adsorption between the dye molecules. Therefore, the lack of surfaces and active sites for adsorption at high concentrations causes less adsorption to be done and reduces the removal efficiency 57 . The simultaneous effect of temperature and initial dye concentration on R% of dyes are seen in Fig. 8e,f. The result indicates the highest removal efficiency took place in higher temperatures with lower dyes concentrations. It can be assumed that www.nature.com/scientificreports/ increasing the temperature increases the kinetic energy of the reacting particles and increases the probability of collisions between dye molecules and the adsorbent surface, as a result, the dye removal efficiency increases 58 . It can be assumed that increasing the temperature increases the kinetic energy of the reacting particles and increases the probability of collisions between dye molecules and the adsorbent surface, as a result, the dye removal efficiency increases, which in line with other previous reports 58,59 .
Optimum conditions of the adsorption process. The desirability function (DF) was used to optimize the values of the independent variables to get the high performance of the ZnS/CuO-CNT in P4R and TA dyes www.nature.com/scientificreports/ removal by applying a scale ranging from 0.0 (undesirable) to 1.0 (very desirable). The profiles for predicted values and desirability function for the maximum R% of dyes are shown in Fig. 9. The dotted red lines indicate the optimum values with a desirability function of 1.0, which were found at solution pH: 2.0, adsorbent dosage: 10 mg, sonication time: 12 min, dyes concentration: 30.0 mg L -1 , and solution temperature: 25 °C. Under the optimum conditions, the predicted R% response for TA and P4R dyes was above 99.0% with a desirability function of 1.0. Also, the experimental R% for TA and P4R were found to be 99.21 ± 2.23, and 98.45 ± 2.54, respectively. Therefore, the predicted results were closely correlated with the experimental data obtained from desirability optimization analysis using CCD-RSM.
Adsorption isotherm models. The adsorption isotherm is critical for describing how the contaminant interacts with the adsorbent and for explaining the adsorbent's adsorption capability. Langmuir, Freundlich, Temkin, and Dubinin-Radushkevich's models were used to analyze the adsorption isotherm. The coefficients of the equations and the isothermal parameters were obtained by calculating the slope and breadth of the regression line in the isotherm diagrams of the adsorption of each dye by the adsorbent from the origin. Table 4 shows the results of these computations as well as the coefficients for the isothermal equations of P4R and TA dyes. According to the results, the Langmuir isotherm is most consistent with the adsorption process for each of the P4R and TA dyes, with a correlation coefficient of 0.999, followed by the Temkin, Freundlich, and Dubinin-Radushkevich isotherms. In the Langmuir model, it is assumed that the adsorbed molecules have no interaction and that the surface adsorption is a monolayer. In the Langmuir relation, K L (L mg -1 ) is the equilibrium constant of surface adsorption, C e (mg L -1 ) is the equilibrium concentration of the adsorbent, q e (mg g -1 ) is the adsorbent's equilibrium adsorption capacity, and q m (mg g -1 ) is the adsorbent's maximum adsorption capacity. The Langmuir isotherm, by defining a dimensionless constant called the separation coefficient (R L ), indicates whether the process of adsorption is desirable or undesirable. For desirable adsorption, the R L value is between 0 and 1. Indicate undesirable surface adsorption, linear adsorption, and irreversible adsorption are all indicated by R L values larger than one, and R L = 1 and R L = 0, respectively 60 . The value of the R L coefficient for adsorption of P4R and TA dyes on ZnS/CuO-CNT were 0.002 and 0.003, respectively. Freundlich isotherm is an isotherm that describes non-ideal and reversible surface adsorption. This experimental model is used for multilayer surface adsorption as well as for heterogeneous surfaces. The K F (L mg -1 ) constant in the Freundlich relation is the absorption isotherm constant, which is related to the adsorption capacity, and n is the Freundlich isotherm constant, which is related to the adsorption severity. If 1/n = 0 the www.nature.com/scientificreports/ adsorption process is irreversible, if 0 < 1/n < 1 the adsorption process is desirable and if 1/n > 1 the adsorption process is undesirable 60 . As can be seen in the table, the value of 1/n for P4R and TA dyes were equal to 0.238 and 0.220, respectively, and indicates the desirability of the adsorption process.
The Temkin isotherm has a factor that includes absorbed-adsorbent interactions. This model considers that the adsorption heat of all molecules on the surface decreases with increasing surface coverage. The quantity of K T (mg L -1 ) in the Temkin isotherm is the constant of the and B 1 (mol J -1 ) is the quantity related to the heat of surface adsorption. The Dubinin-Radushkevich isotherm is an experimental model. In the Dubinin-Radushkovich model, if the amount of adsorption energy (Kj mol -1 ) is less than 8, the adsorption is of the physical type due to weak van der Waals forces, If it is in the range of 8-16, the adsorption follows the ion exchange mechanism, and If E is between 20 and 40, chemical adsorption has taken place 61,62 . The E value calculated in the Dubinin-Radushkovich isotherm for TA dye was 9.054 kJ mol -1 and shows the presence of ion exchange mechanism (E > 8 kJ mol -1 ). For the P4R dye, the E value was 6.704 kJ mol -1 , suggesting that P4R dye has been adsorbed mainly via physisorption (E < 8 kJ mol -1 ).
The nonlinear Chi-square (χ 2 ) test was used to assess the fitness of isotherm models in addition to the correlation coefficient of determination (R 2 ). If the data from a model and the experimental data are comparable, χ 2 is a small number; if they are significantly different, χ 2 is a large number 63 . As can be seen, the Langmuir model exhibited lower χ 2 values than other isotherms, which lead to a better fit of the calculated data from the Langmuir model and experimental results.
Adsorption kinetic studies. The various kinetic models were used to study the adsorption kinetics (Table 5). According to the First-order-kinetic model, infiltration occurs inside a boundary layer. This model is based on the adsorbent capacity, which states that variations in the quantity of adsorption over time are proportionate to the number of vacant places on the adsorbent level. The pseudo-second-order-kinetic model is based on solid-phase adsorption and indicates that chemical adsorption is the primary control mechanism in the adsorption function. It also places that the chemical adsorption steps slow down the adsorption process. Weber-Morris intraparticle diffusion model is used to study the mechanism of the adsorption kinetics. The quantitative adsorption process may be controlled by one or more steps, including external penetration or penetration into the film, penetration into the pores, penetration into the level, and adsorption on the surface of the pores, or a combination thereof. In a solution that is stirred rapidly, mass transfer through diffusion can be expressed by an apparent diffusion coefficient calculated from experimental data 64,65 . Elovich model is also one of the most widely used models and unlike previous models that are compatible with physical adsorption, it shows the chemical interaction between adsorbate and adsorbent 66 . According to the given models, the factors calculated from each model are reported in Table 5. In these equations, q e (mg g -1 ) and q t (mg g -1 ) are the adsorption capacities at equilibrium time and time of the t, respectively. k 1 (min -1 ), k 2 (g mg -1 min -1 ), and h (mg g -1 min -1 ) are the first-order velocity constant, second-order velocity constant, and initial sorption rate, respectively. K diff (mg g -1 min 0.5 ) is the penetration rate constant, and C (g mg -1 ) is a constant that gives an idea of the thickness of the boundary layer. Also, α and β are the indicators of the Elovich model 67 . The results show that the adsorption process follows the Pseudo-second-order-kinetic model based on the R 2 value of P4R (0.996) and TA (0.996) dyes. Also, the amount of equilibrium capacity calculated (q e, calc ) from the Pseudo-second-order-kinetic model is more consistent with the equilibrium capacity obtained from experiments (q e, exp ). Therefore, a better fit www.nature.com/scientificreports/ of the results with second-order kinetics model indicated that the chemical adsorption was the primary control mechanism in the adsorption process.
Although the adsorption process of TA and P4R by ZnS/CuO-CNT was followed well by the Pseudo-secondorder-kinetic model, however, the diffusion mechanism and the rate-controlling step in the adsorption process are not still confirmed. Therefore, the Weber-Morris model was used to study the actual rate-controlling steps in the adsorption process. The Weber-Morris fitted plot of qt against t 1\2 is shown in Fig. 10. The result illustrates a multilinear plot in three sections, which confirmed that the adsorption of TA and P4R by ZnS/CuO-CNT was controlled by a three-stage mechanism. The first stage of adsorption had a higher slope and up to 1.8 (min) 1/2 , which was related to the boundary layer diffusion of dye molecules onto the surface of the adsorbent. Then, the slope of the straight section decreased gradually in the second stage and showed a gradual linear adsorption fit www.nature.com/scientificreports/ between 1.9 and 2.8 (min) 1/2 , which was related to the diffusion of dye molecules into the interior pore structure of the adsorbent. Finally, the slope of the straight section approached almost zero in the last stage of adsorption from 2.9 to 3.2 (min) 1/2 , which confirmed that the adsorption process reached the equilibrium condition. In this stage, all active site of the adsorbent was wholly saturated with dye molecules. Also, the non-zero values of Ci at all the stages showed that the adsorption process was not only intra-particle diffusion, and other interaction mechanisms including adsorption on the external surface might be acting simultaneously to it 68 .
Thermodynamic parameters. In the adsorption process, the determination of thermodynamic parameters is important to detecting the endothermic or exothermic of the adsorption process as well as for determining its spontaneity. In thermodynamic studies of the adsorption process, it is necessary to determine three parameters, standard enthalpy (∆H°), standard Gibbs free energy (∆G°), and standard entropy (∆S°) 67 . Thermodynamic variables at different temperatures were investigated and the results are collected in Table 6. The positive values of ∆H°, indicating that the adsorption process of the dyes by ZnS/CuO-CNT was an endothermic reaction. Also, positive values of ∆S° indicate an increase in irregularity in the solid-soluble surface during the adsorption process. In other words, the positive entropy changes of the standard system indicate an increase in irregularities in ZnS/CuO-CNT in the adsorption process of P4R and TA dyes compared to the initial state before the adsorption process. Negative changes in standard Gibbs free energy showed that the adsorption process of P4R and TA dyes was spontaneous. Also, it has become more negative with increasing temperature, which indicates that the spontaneity of the reaction has increased at higher temperatures.

Proposed adsorption mechanism.
To better understand the possible adsorption mechanism, isotherm data, and adsorption kinetics were examined and showed that the adsorption of P4R and TA dyes by ZnS/ CuO-CNT can occur due to both physical and chemical reactions as follows: (i) the slope (1/n) of Freundlich isotherm was less than one for both dyes and indicated the main process of the adsorption was chemisorption.
(ii) The adsorption energy (E) of Dubinin-Radushkevich isotherm, indicating that the adsorption process of TA and P4R dyes exhibit chemical and physical adsorption behaviors, respectively. (iii) better agreement of the adsorption data with the Pseudo-second-order-kinetic model shows that the chemisorption is the dominant control mechanism for both dyes. Also, the FTIR spectra of ZnS/CuO-CNT before and after adsorption of TA and P4R dyes are present in Fig. 11. The result showed that the C-O, C = C, − COOH, C-H, O-H, CuO, and ZnS functional groups of ZnS/CuO-CNT shifted after the adsorption process, which could indicate the involvement of hydrogen bonding, π-π, n-π, and electrostatic interaction mechanisms in the adsorption process. On the other hand, the appearance of absorption peaks related to the TA and P4R dyes in the after adsorption FTIR spectrum of ZnS/CuO-CNT also indicates the presence of chemisorption process 69 .

Adsorption of TA and P4R food dyes in real effluents.
To evaluate the adsorbent performance in real conditions, the simultaneous removal of P4R and TA dyes was investigated using ZnS/CuO-CNT adsorbent in distilled water and also real effluents including tap water, river water, and polluted water. Figure 12 shows the characteristics of the effluent following adsorption with the ZnS/CuO-CNT adsorbent. As can be observed, adsorbent has good effectiveness in removing binary combination dyes at the same time and was able to remove over 90% of dye residuals. The results demonstrated the ZnO/CuO-CNT can be very successful in treating dye effluents in real environmental conditions. Adsorbent regeneration and reuse experiments. The regenerate ability and reused the spent adsorbent is very critical in economical applications. For this purpose, the dye-loaded adsorbent was regenerated by using the desorbing agent of NaOH (0.1 mol L −1 )/ethanol at ambient temperature, washed with ultrapure water, and after being dried in the oven (120 ºC for 1 h) was reused in the next cycle of adsorption experiments. Figure 13 shows the experimental results of the dyes removal efficiency by the ZnS/CuO-CNT adsorbent in each adsorption-desorption cycle. As can be seen, absorbent reuse was performed in six stages of the adsorption/ desorption cycles. The results showed that the regenerated adsorbent after the five adsorption-desorption recycles could still remove dyes more than 90% from the aqueous solution. Further, the relatively high desorption efficiency and good regeneration of adsorbent by desorption agent indicate that the adsorptive removal of dyes www.nature.com/scientificreports/ by the adsorbent may have mostly occurred by electrostatic interactions 70 , which corroborates well the explanation of the initial solution pH effect on dyes removal efficiency investigated in section 'Simultaneous effect of independent variables on dyes removal efficiency' .

Comparison of maximum removal capacity with different adsorbents.
One of the things that should be considered in the performance evaluation of an adsorbent is the study of the maximum adsorption capacity (Q max ) of the prepared adsorbent compared to other common adsorbents. Table 7 shows the results of the comparison Q max of the ZnS/CuO-CNT adsorbent with other adsorbents. As can be seen, the ZnS/CuO-CNT has a good Qmax for both P4R and TA dyes and is competitive with other reported adsorbents in terms of performance.

Conclusion
In this study, the ZnS/CuO-CNT nanocomposite was synthesized by the in-situ hydrothermal synthesis, and its performance as a potential adsorbent in the simultaneous ultrasound-assisted adsorptive removal of a binary mixture of P4R and TA dyes from contaminated waters was systematically investigated. The physico-chemical characteristics of synthesized adsorbent were studied by XRD, FESEM, FTIR, BET, and ZP analysis. The  www.nature.com/scientificreports/ independent process variables including adsorbent dosage, initial solution pH, sonication time, initial dyes concentration, and temperature were considered and statistically optimized to maximize P4R and TA dyes adsorptive removal efficiency using the CCD-RSM design. Based on ANOVA analysis, the experiment was systematically adopted with the confidence level of p < 0.05, and the statistically significant model was developed for dyes removal through regression analysis (R 2 = 0.99). The optimal values were adsorbent doseage of 10 mg, pH value of 2.0, sonication time of 12 min, dyes concentration of 30.0 mg L −1 , and solution temperature of 25 °C. The results showed that the removal efficiency increased at acidic pH and is also directly related to the sonication time, adsorbent dosage, and temperature. In contrast, increasing the concentration of dyes reduces the removal efficiency. Under the optimum conditions, the predicted R% response for TA and P4R dyes was above 99.0% with a desirability function of 1.0. Also, the experimental R% for TA and P4R were found to be 99.21 ± 2.23, and 98.45 ± 2.54, respectively. The optimum conditions were applied to study the adsorption isotherms, kinetics and thermodynamics. In addition, the adsorbent reusability results revealed that the adsorbent after the five adsorption-desorption re-cycles could still remove dyes more than 90% from the aqueous solution.  www.nature.com/scientificreports/